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1. Introduction 

Periodic nanostructured media, also known as photonic crystals (see, e.g., [1,2] and 
references therein), have put forth a wide range of possible applications in optical 
communication, optoelectronics and optical means of data transmission. In addition, 
photonic crystals have inspired a lot of fundamental research on light-matter interaction, 
since electromagnetic wave propagation phenomena in strongly inhomogeneous media 
have been for the first time available for direct theoretical investigation. 

Moreover, during the recent decade it was shown that further alteration of topology 
of inhomogeneous media leads to an even wider choice of optical materials. Within 
this scope, non-periodic deterministic (NPD) structures have been extensively studied 
during the recent decade. Namely, quasiperiodic [3-5] and fractal [6-8] nanostructures 
are known to possess distinct optical properties, not present in either periodic or 
disordered systems. The most interesting among these are spectral self-similarity and 
critical eigenstates in Fibonacci quasiperiodic lattices [3] and spectral scalability in 
Cantor fractal multilayers [6], which was shown to directly result from geometrical self- 
similarity [8]. 

However, nearly in all the research in this field only optically isotropic constituent 
materials are considered. On the other hand, there have been recent reports with 
investigation of anisotropic materials embedded into a periodic structure. For example, 
it was shown that anisotropy can cause even a one- dimensional (ID) multilayer structure 
to exhibit negative refraction [9, 10], a phenomenon earlier attributed exclusively 
to higher- dimensional systems, as well as omnidirectional total transmission [9, 11]. 
In addition, interesting results have been recently obtained with weakly modulated 
periodic structures made of anisotropic and gyrotropic materials [12]. Extensive ongoing 
analytical research of anisotropic and gyrotropic multilayers [13-15], as well as studies of 
a simpler case of chiral media in 2D [16] and 3D [17] photonic crystals, is also restricted 
to periodic systems. 

From what has been said, there is a need to "bridge the gap" and combine anisotropy 
and non-periodicity in the same structure. It is of particular interest to find out how 
optical anisotropy of constituent materials can affect the spectral properties caused by 
deterministic non-periodic geometry. 

Our earlier communication [18] showed that for any multilayer structure exhibiting 
sharp resonance peaks inside a band gap (generally, a characteristic feature of NPD 
geometry) anisotropy results in polarization- induced peak splitting (figure [T^i,b). In this 
regime, the structure acts as a nanosized absorptionless frequency-selective polarizer, 
which can in itself find certain applications. Moreover, in the case when a doublet 
undergoes polarization-induced splitting in such a way that the components are matched 
(figure[lfc), the structure acts as a quarter- wave retarder (QWR) and converts a linearly- 
polarized beam into circularly-polarized and vice versa. The size of such a retarder is 
smaller than that of a bulk QWR, and the content of birefringent material is more 
than 10 times decreased, which suggests enhancement of effective birefringence when 
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anisotropic material in organized in deterministic non-periodic manner. 

Of the geometries in question, fractal structures are of special interest because their 
spectra can provide resonance peaks with controlled multiplicity due to sequential peak 
splitting [6]. 

In this paper, we move on to further analyze the optical properties of non-periodic 
media containing optically anisotropic constituent materials. Along with conventional 
birefringence, we analyze the effects caused by optical activity (or gyrotropy) [19]. We 
show numerically as well as analytically that in the case of bi-isotropic multilayers 
(isotropic optical activity introduced into an initially isotropic medium) gyrotropy does 
not cause polarization-induced splitting - contrary to what might have been expected. 
In fact, gyrotropy even does not manifest itself in the spectral characteristics of the 
multilayer, its sole influence being a geometry-independent polarization plane rotation 
of the transmitted wave. Furthermore, we show that optical activity can be made use 
of to modify the polarization sensitivity and the output polarization of a birefringent 
polarizer described above, meaning both polarization orientation and ellipticity. Possible 
applications for integrated optics and optical engineering are also outlined. 

To account for gyrotropy, we make use of coordinate-free operator formalism for 
the general case of bi-anisotropic multilayers [20,21]. It is outlined briefly in Sectional 
The simplest case of a multilayer with bi-isotropic layers is analytically and numerically 
investigated in Section El It is shown that gyrotropy can only provide polarization plane 
rotation completely independent of structure geometry. The combination of isotropic 
optical activity and uniaxial birefringence is analyzed in Section HJ and it is shown that 
in this case gyrotropy can be used to control the output polarization of an integrated 
optical polarizer. Section El summarizes the paper. 
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Figure 1. A portion of spectrum of a 27-layer Cantor fractal multilayer structure as 
a typical spectrum of an NPD structure (a) along with the image of one peak (b) 
at different orientations of input wave polarization. The peak exhibits polarization- 
induced splitting [18]. Sometimes, a doublet may split in such a way that matching 
occurs (c). The inset shows the structure geometry, which can be written as 
gagaaagagaaaaaaaaagagaaagag , g being GaN {e — 5.685, e e — 5.892) and a being air 
quarter- wave layers. Frequency is normalized to luq defined ELS Tlgdg — Tlada = A /4 = 
ttc/2uj [18]. 
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Figure 2. Illustration of an ambient medium / layer interface: on the definition of a 
natural basis in J2J). 

2. Evolution operator for anisotropic and gyrotropic media 

Before proceeding, it is necessary to briefly touch upon the theoretical approaches used. 
Contrary to the case with isotropic media when common transfer matrix methods 
with scalar field values would suffice, anisotropic and gyrotropic media cannot be 
regarded without accounting for polarization effects. Such effects include, for instance, 
polarization coupling at layer interfaces, so in the general case wave polarization cannot 
be separated into independent states and full vector calculations have to be performed. 

Among the methods for such calculation, we have used the coordinate-free operator 
method (also known as Fedorov's covariant approach [19,20,22]), developed for stratified 
systems in [21]. In this method most equations remain in compact form containing 
only the inherent parameters of the structure. This provides flexibility and facilitates 
analytical investigation of the equations used. 

To account for arbitrary anisotropy and gyrotropy, the following form of material 
equations is used [21]: 

D = eE + zaH, B = /iH + i(3E. (1) 

Here e and /i are dielectric permittivity and magnetic permeability, respectively. In 
optically anisotropic media, e is a tensor rather than just a scalar coefficient. In materials 
possessing magnetic anisotropy, \i can also be a tensor. The quantities a and j3, also 
tensorial in the general called gyration pseudotensors and are responsible for 

the material's optical activity. In non-absorbing media they must satisfy the relation 
(3 = —ofi (Hermitian conjugation). The material equations ((TJ can in theory encompass 
all possible cases of anisotropy and gyrotropy [19]. 

Considering a plane wave incident from isotropic ambient medium with refractive 
index n onto a slab of a material with parameters given in ([TJ, one can define the 
wave by its vector of refraction m = uqii, n being the vector of wave normal (figure |2j). 
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Analogously, the boundary can be characterized by the unit normal vector q and the 
corresponding projection operator 



q ® q 



(q 



x\2 



which projects any vector into the interface plane. Here the notation u <g> v denotes 
an outer product (also termed the dyad [21]) between two vectors u and v, which is 
a tensor with elements in the form (u ® v)y = UiVj, and the symbol q x stands for an 
antisymmetric tensor defined via vector cross product as (q x )u = q x u, or via fully 
antisymmetric Levi-Civita's pseudotensor Syk as (q x )j k = £ijkQi- 

Using these elements as a reference, one can further introduce the vectors 

b = im, a = b x q. (2) 



The three vectors (ao = a/|a|,bo = b/|b|,q) are of unit length and mutually 
orthogonal. Thus they form a basis in space (see figure EJ). Note that at normal incidence 
a = b = 0, which results from axial symmetry of such a system. In this particular case 
only one of the basis vectors can be introduced naturally, while a and b are subject 
to arbitrary choosing. 

It was shown earlier [21] that Maxwell's equations with Q can be transformed into 
an operator differential equation 



d W{x) 



dx 



ikMW{x) 



W(x) 



IH(x) 
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, M = 
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D 



(3) 



Here M is a characteristic block matrix (or characteristic operator) whose elements are 
operators dependent on the material parameters of the slab. It can be represented by 
a 6 x 6 matrix. However, the block matrix notation allows to discriminate between real 
space [where we have introduced the vectors (J2J] and "artificial" 2x2 space related 
to block vectors W, introduced because electric and magnetic field evolution has to be 
accounted for simultaneously. To differentiate from real-space vectors and operators, 
block vectors and matrices are denoted by light capital symbols. 

So the components of the characteristic matrix in Q are operators acting in real 
space. For their full form one can refer to [21]. Throughout this paper, normal incidence 
will be considered both in analytical and in numerical investigations. This makes certain 
sense from the experimental point of view because otherwise gyrotropy is known to be 
completely shadowed by common birefringence. In this case, 



A = iq x al + [q x eq <8> v 3 + iq x aq (g> Vi] , 

B = -q x eq x + [q x eq ® (q x v 4 ) + iq x aq <g> (q x v 2 )], 

C = IfiI + [H(3q <g> v 3 + Ifiq <g> vi], 

D = -Uf3q x + [Upq ® (q x v 4 ) + i>q ® (q x v 2 )], 



(4) 
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Here the four auxiliary vectors v& are introduced as follows. 

vi = -Qq(e q /i + fi q a)I, v 2 = iQq((3 g e - e q (3)I, 



(5) 



V3 = iQq(a q [i - fi q a)I, v 4 = —Qq(ji q e + a q (3)L 

In these expressions, the subscript q means multiplication by q on both sides 
(s q = qeq, etc.), Q = (s q fi q + being the normalizing factor. 

Having established the characteristic operator (J3J for a layer with thickness d, one 
can proceed to solve the equation Q by taking the matrix exponential 



W(0)=PW(d), P = exp 



(6) 



Here the operator P is termed the propagator or evolution operator. It completely 
determines the wave propagation in one layer made of any material. Once known for 
each layer of an iV-layer structure, the propagator for the whole system can be obtained 
through operator multiplication: 

P = P 7V P 7V -i---Pi. (7) 

The resulting evolution operator P can be used to establish the reflection, 
transmission, and spatial field pattern for the structure in question. The details are 
given in [21]. 

Note that the form of (jUJE)) is quite general, and can be simplified greatly in some 
cases. For instance, for non-gyrotropic isotropic media (at = (3 — 0, e and \i are scalar) 
after some algebra one can reduce (jlj) to 



T(0 



el 
til 



(8) 



Note that we have used the symbol e rather than e to distinguish between scalar 
and tensor permittivity, respectively. To avoid confusion, this will be done throughout 
the article. In what follows, the rest quantities (/x, at, (3) will be assumed scalar unless 
explicitly stated otherwise. 

As a final note, one can see from (J3J) that the block vector W contains only in-plane 
components of the field vectors. Therefore in the case of normal incidence and isotropic 
ambient medium with refractive index n 



(9) 



m 
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q x E 




7H 



7 



1 

n 



a ® a + b ® b c 



-I. 



n 



where 7 is called a surface impedance tensor and can be derived for this case from [21]. 
Note that once we consider finite-sized multilayers, 7 depends only on the properties of 
ambient medium and thus always has the form since the resulting propagator 
acts in (0) at the points just outside the structure boundaries. 
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3. Multilayer structures with bi-isotropic layers 

In this section, we consider the simplest case of gyrotropic materials, when all the four 
tensors (e, fi, a, (3) are scalar. Such parameters correspond, for instance, to cubic 
crystals or optically active liquids, which are sometimes called bi-isotropic or chiral 
media. Because of their relative simplicity, these media were subject to earlier general 
analytical studies both as bulk media [22, 23] and as stratified media or multilayers 
[24-32]. However, as noted above, the latter was either focused on general problems of 
a chiral layer interface [24-27] or concerned with periodic structures [24,28-31]. Even in 
the recent, more general theoretical works [32] the authors still confine their numerical 
studies to periodic media with bi-isotropic layers. 

Similarly to uniaxially birefringent ones, bi-isotropic media are known to change the 
effective refractive index depending on the beam polarization, but with circular rather 
than linear eigenpolarizations (see [27] and references therein). Accordingly, one may 
expect that an NPD multilayer with bi-isotropic layers will exhibit similar polarization- 
induced splitting [18] and act as a polarizer for circularly polarized light. 

Numerical experiments, however, show that this is entirely not the case. We have 
used the same system as in [18], a 27-layer GaN/air fractal structure shown in the 
inset of figure All GaN layers are artificially rendered isotropic by setting e D = e e . 
Due to absence of birefringence, polarization-induced splitting vanishes, so only one 
transmission peak remains and the spectrum becomes that for the extraordinary wave 
in figure EJo. Variable isotropic optical activity is then introduced by imposing scalar 
non-zero a (ranging up to 0.2, which is somewhat unrealistically high but chosen so as 
to fully cover the range currently available in experiments), so all GaN layers become 
bi-isotropic. 

As a increases, the transmission spectrum shows no modification whatsoever. The 
reflected wave does not change its polarization, either. The transmitted wave, however, 
exhibits polarization rotation. Figure El shows polarization diagrams for the transmitted 
wave for different input polarizations and for different values of a, for two incident 
wave frequencies. These polarization diagrams is what a photodetector would register if 
placed behind a rotating analyzer in front of the output beam. It can be seen that this 
rotation is the same both at the transmission peak (figure OH-f) and off-peak (figure EK- 
c). So the polarization rotation is uniform, i.e., it has a similar manner at all frequencies 
and regardless of the spectral features. The angle of rotation can be written as 



So it can be concluded that the influence of gyrotropy in this case is completely 
incoherent, i.e., totally independent of the geometrical properties of the multilayer. 

To illustrate this seemingly unexpected result, let us recall that in bi-isotropic 
media the polarization rotation direction is known to be reversed after normal-incidence 
reflection (so-called cross-polarized case described in [26]), in full accordance with the 



N 




(10) 
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fact that the medium itself does not exhibit any specific directions. So the reflected part 
of the incident beam (see figure will come to the initial interface (0|1) at its incident 
polarization, and this is true regardless of the number of multiple reflections taken into 
consideration. Similarly, at the second layer interface (1|2) all beams will have the same 
polarization, rotated with respect to the incident wave polarization according to (|10|) . 
Since the transmittance is determined by interference effects resulting from multiple 
reflections, optical activity in a bi-isotropic layer cannot make any contribution except 
the above-mentioned polarization rotation (figure EJ). 

This explanation should hold regardless of the number of layers in a multilayer 
structure or of the structure geometry, because interference of multiply reflected beams 
should occur likewise at all the layer interfaces provided all the constituent layers are 
either isotropic or bi-isotropic. This was shown for periodic stacks in [29] and can also 
be derived from single-layer and single- interface results obtained in [24,26]. 

To provide a consistent back-up for arbitrary NPD multilayers, let us calculate the 
evolution operator for a bi-isotropic layer. Starting with (j^EJ), one can immediately 
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Figure 3. Polar plots of transmitted wave and incident wave polarization diagrams 
(namely, the dependence of |(x ® x) H| (tp) where the vector x(^) = a cos ip + b sin i\j 
rotates in the plane of interface) for a 27-layer fractal multilayer containing bi-isotropic 
layers with variable optical activity a. The plots are given for off-peak [w/wo — 1.0930, 
(a-c)] and on-peak frequencies [lu/ujq = 1.0938, (d-f)]. The incident wave polarization 
varies between 0° (a), 22.5° (b), 45° (c,f), 67.5° (e), and 90° (d). 
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Figure 4. The illustration of spectrally uniform polarization rotation in bi-isotropic 
media. First, the incident wave hits the left interface (a). As the wave passes through 
the layer, its polarization is rotated (b). At the right interface (c) the wave is partially 
transmitted and partially back-reflected. As the back-reflected wave travels back, its 
polarization is rotated back (b), so it arrives at the left interface at exactly the same 
polarization as the incident wave (a). The process is repeated as this wave is partially 
reflected and partially transmitted again. (Double arrows denote propagation direction 
of the waves shown by the direction of their H- vectors). 



notice that for all the auxiliary vectors (jSJ) the expressions in brackets are scalars and 
thus commute with q, leaving the product q/, which is identically zero. Remembering 



that Jq x = 
reduced to 



q x J = q x and that the absence of absorption implies f3 



-a, 



M (b) 



el 
zaq > 



can be 



Taking the matrix exponential, one can, after some algebra, find the propagator to 
equal \ 



cos(kdy/eJI)P(a 



sin(kdy/eJI)P(a 



sm(kdy/ejI)P(a) cos{kdJep)P{a 



(12) 



P(a) = [cos(kda)I — sin(£;<ia)q > 



where k = u/c. 

The evolution operator for the isotropic layer with the same e, /x, and d can be 
found either by imposing a = or directly from (j8j), which results in 



p(0 



cos(kdy/ejI)I -j== sm(kdyfejl)I 
cos(kdyfql)I 



^=sm(kd^jl)l 



(13) 



X If we are working in a 3D real space, the propagator will formally contain an additional term q <£> q 
added to diagonal terms. However, in our system only normal incidence is considered. Hence for all 
field vectors qH = 0. Since, naturally, 7q — = q x q, all calculations are unchanged if performed in 
a 2D subspace associated with the interface plane. Formally speaking, we can consider blocks of M to 
be 2 x 2 rather than 3x3 matrices before taking the exponential. 
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We can see that every real-space block of both ()T2|) and (jl3j) has identical 
geometrical structure, with different coefficients for different components of the block 
vectors W in (El). This can be written as 



P (fc) = (cos(kda)I - sin(kda)q x ) L, P (i) = IJL 



where 



L 



cos(kd^/eJI) -^=sm(kdy/ejl) 



sm(kd^yejl) cos(kdy/eJI) 



(14) 



(15) 



Here one can note that the block matrix L consists of scalar quantities, so in terms 
of real space it is a coefficient matrix, which does not affect polarization properties. This 
essentially results from normal incidence. Also note that this coefficient matrix is the 
only part in P that depends on propagation phase ip = ^/epduo / c. 

Now let us show what such a propagator does to a normally incident vector H . 
For simplicity and without loss of generality, we can assume H = H sl , after which it 
can be seen from (J^J) that propagation in isotropic layer results in 



W = P W W = H Ih 
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—an 

no u 



HnL 



1 

no 



Jan = Hn 



t 
f 



a . 



(16) 



One can see that there is no polarization rotation. On the other hand, substituting 
the evolution operator of a bi-isotropic layer analogously yields 



W = P (b) W n = Hn 



t 
f 



(ao cos(kda) — b sm(kda)) 



;i7) 



We can see that the polarization of the transmitted wave is rotated by the angle 
of (ft = kda, while the values of t and t! (which contribute to the transmittance) remain 
the same. Besides, it can be stated that the propagator structure for the bi-isotropic 
layer in (fPfj) explicitly contains the value of 0. 

To advance from a single layer to a multilayer structure, let us investigate what 
happens when propagators in the form ((Hj) are multiplied. Naturally, we assume 
different e, fi, d, and a (if any) for the two adjacent layers. 

So, a product of two isotropic propagators reads 



P 



12 



(LxL 2 ) / = L 12 J 



(18) 



(Li/) (L 2 J) 

One can see that while the coefficient matrices Lj are multiplied (and hence the 
complex spectral effects are seen in isotropic multilayer structures), the real-space 
structure of the propagators remains proportional to J, i.e., polarization-independent. 

Likewise, a product between isotropic and bi-isotropic propagators reads 



P 



iib) 



12 



pWpW 



L 12 (cos(fc<i 2 Q; 2 )/ — sin(/c<i 2 a 2 )q x ) (19) 

Again, the coefficient matrices are multiplied in the same manner, while real-space 
geometrical structure is preserved, so the angle of polarization rotation is unchanged 
compared to that for a single bi-isotropic layer. 
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Finally, taking a product between two bi-isotropic propagators, after some algebra 
we arrive at 

— (cos{kdiot.i + kd 2 a 2 )I — sin(A;(i 1 a 1 + kd 2 a 2 )q x ) ■ (20) 

So we see that the coefficient matrices are again multiplied likewise, and hence 
no change to the spectral properties provided that the same Li and L 2 are used 
in ()18|) and (|20|) - which is true if the corresponding layers have the same e, fx, and d. 
Looking at the geometrical structure of the propagators, we see that 

012 = 0i + 02 = — (aidi + a 2 d 2 ) . (21) 
c 

Certainly, as the propagator for any multilayer is built according to (J7j). and any 
multiplication according to (p"8H20|) leads to the fact that spectra are independent of a 
while polarization rotation sums up, so (|2ip is directly generalized to (|10|). Therefore 
the reasoning applicable to a single bi-isotropic layer holds for any fractal structure (or, 
for that matter, for an arbitrarily designed binary multilayer composed of isotropic and 
bi-isotropic layers). 

So far we have considered a linearly polarized incident wave. If the polarization is 
circular, it can be shown that change in effective refractive index does occur, however 
it does not manifest itself in the transmission spectra. Detailed explanation is provided 
Appendix A 



in 



To summarize, we have shown analytically that in multilayers consisting of isotropic 
and bi-isotropic layers in arbitrary combination the optical spectra are exactly the same 
as if there were no optical activity whatsoever, regardless of incident wave polarization. 
Instead, optical activity rotates the transmitted wave polarization, and the amount of 
this rotation is totally independent of the multilayer's geometrical composition and is 
described by a simple equation ([TT)]) . 

This result can also be understood from the point of view that in a bi-isotropic 
multilayer system the multiple-reflected beams will interfere at each layer interface in 
exactly the same way as they do in the isotropic multilayer with the same e, li, and d. 
Since the optical spectra are determined by the nature of such interference, they are 
unmodified despite polarization rotation caused by each bi-isotropic layer and resulting 
in overall uniform polarization rotation. 

Note that the analytical derivation in this case is rigorous and exact, without any 
assumptions on the values of parameters used. The only assumption regards the general 
applicability of local material equations (0), which restricts the layer thickness to the 
wavelength such that d ^> 0.001A [19] - a condition well met in the current state-of-the- 
art multilayers in the optical range. 

In conclusion to this section, let us point out that the uniform and incoherent 
nature of the polarization plane rotation allows to use bi-isotropic materials to rotate 
polarization in devices with any geometry, without risk of additional interference 
effects. It can also be possible to combine spectral and polarization-rotating device 
in one multilayer, which should be beneficial for miniaturization of integrated optical 
components. 
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Figure 5. The dependence of transmittance on frequency and incident wave 
polarization for a 27-layer GaN/air fractal structure {inset) without (a) and with 
gyrotropy (b) with a = 0.2. Lighter areas denote higher transmittance. 



4. Gyrotropy and polarization-induced splitting 

Gradually complicating our structure, let us now consider a birefringent medium with 
a uniaxial permittivity tensor 

e = t + (e e - t ) c ® c (22) 

where, as noted before, e and e are used to discriminate between tensor and 
scalar permittivity, respectively, and subscripts o and e stand for "ordinary" and 
"extraordinary". A unit vector c determines the optical axis orientation. In multilayer 
structures under study, all uniaxial layers have similar orientation of c such that cq = 0. 
Without loss of generality we can assume c = bo- As before, we assume no magnetic 
anisotropy, and isotropic optical activity. 

Substituting (J22|) into (JIHHJ) lets us arrive after some simple algebra at 



M (a) 



iaq* el + iQfi q q x eq (g> q x qel 
\il iaq x 



iaq 



el 



fil iaq y 



Without optical activity (a 
propagator can be found to equal 



0) and within our assumptions (c 
P (a) = L o a ® a + L e b g> b 



(23) 
bo) the 

(24) 



where 



sm(kdy/€jn) 



■ sm(kd^/€jfi) 



cos(kd^/€jfi) 



o, e. 



(25) 



It can be easily seen that such a layer is able to exhibit polarization-induced peak 
splitting in the spectrum. Indeed, when a field vector is oriented along a (or, more 
generally speaking, perpendicular to c), it will only interact with one term in the 
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propagator, namely, the one containing L , and will propagate exactly as if all uniaxial 
layers were isotropic with e = e D [compare (fT3|) and ([25)1] . Similarly, when a field 
vector is oriented along bo (parallel to c), the same would be true for L e , and again 
all uniaxial layers will effectively be isotropic, this time with e = e e . So, one can see 
that if the incident wave is polarized along either of the eigenvectors of P (it is then 
called an eigenwave), it will propagate in exactly the same way as if all uniaxial layers 
were isotropic with some effective refractive index. Hence the spectra will be similar in 
shape, the only difference associated with a change in effective e, which in turn causes 
all spectral features to shift in frequency. 

Unfortunately, in presence of optical activity the tensorial nature of e results in a 
very complicated structure of the propagator, which, while capable of being evaluated 
symbolically, is nearly useless for further analysis. 

It can be shown, however, that if we assume both a and 5e = \ |e e — e | to be 
small, an approximate model propagator can be used. The details on its derivation 
and applicability are covered in |Appendix B| Reduced with respect to 5e with the 



product a5e neglected, it reads 

P (9) « L (cos{kda)I - sin(fcda)q x ) + 5eh' (a <g> a - b ® b ) . (26) 

Here the coefficient matrix L is the same as that for isotropic media [compare 
()B.9|) and (JT3J)] if we substitute e = \ (e e + e D ), i.e., if an average between ordinary 
and extraordinary dielectric constant is used as an effective value of e. The coefficient 
matrix L' is explicitly written as (jB.10|) . 

On the other hand, one can reduce the same model propagator with respect to a, 
in which case 

a L o a (g> a + L e b <g> b - kdahq* , (27) 



where the matrices L G and L e are the same as those without gyrotropy [see (|24H25j) ]. 
and the additional coefficient matrix L given in (jB.9|) is the same as for isotropic case 
in JTSJ). 

Looking at (}26|) and (|2*7jl . one can see that even within the bounds of the 
approximation used, the behavior of the multilayer in question changes dramatically. 
Compared to a bi-isotropic evolution operator (^3J), there is a symmetric addition 
proportional to 5e, which causes perturbation to polarization plane rotation, making 
the propagator polarization-sensitive. Hence, the propagator multiplication is no longer 
described by (J20|) . and the angle of rotation becomes dependent on the incident wave 
polarization. 

Compared to a birefringent evolution operator (}24|) . there is an antisymmetric 
addition proportional to a (and, more generally, this a comes from decomposition 
of sin kda). We know from (|T4"|) that such a term in the propagator is responsible 
for polarization plane rotation. 

As it is not easy to decide in favor of one on these two interpretations, it is safe 
to assume that one layer made of material in question will exhibit polarization rotation 



Spectral and polarization effects in deterministically nonperiodic multilayers containing optically anisotropic 

combined (more or less on equal terms) with spectral properties of birefringent multilayer 
structures. 

Numerical results show that this is indeed the case. For simulation, the same 
GaN/air fractal multilayer as in the previous section was used (shown again in figure EJ), 
but without rendering the GaN layers isotropic. Optical activity was likewise artificially 
introduced into GaN, a ranging up to 0.2. 

As one can see in figure there is no change in the nature of polarization- induced 
splitting, and the frequencies of the doublet components are not changed, either. This 
agrees with ()27|) where we see the coefficient matrices L Q and L e to be the same regardless 
of optical activity. However, the sensitivity of splitting to input polarization has been 
modified dramatically. Without gyrotropy there is a clear symmetric picture with peak 
maxima corresponding to the cases when the incident wave coincides with one of the 
eigenpolarizations (figure EH). When gyrotropy is present, the peak maxima become 
shifted by about 30 degrees, and the transmittance at maximum does not reach unity 
(figure Eb)- 
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Figure 6. Polar plots of transmitted wave and incident wave polarization diagrams 
(the same as in Fig. [3| for a 27-layer fractal multilayer containing birefringent layers 
with variable isotropic optical activity a. The plots are given for the same incident wave 
polarizations for both polarization-split peaks [oj/ujq = 1.0923 (a-c), oj/ujq = 1.0938 
(d-f)]. The incident wave polarization varies between 0° (a), 22.5° (b), 45° (c,f), 
67.5° (e), and 90° (d). 
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Figure 7. The dependence of polarization rotation angle 4> on the value of a for 
isotropic (a) and birefringent (b) multilayer structure. The incident wave polarization 
is according to Fig. Et for the split peak frequencies uj/oj = 1.0923 (dashed) and 
uj/ujo = 1.0938 (solid). The inset shows a slight difference in the slope for the two lines 
according to (fTTTf) - 



To proceed with in-depth analysis, let us look at the transmitted wave polarization 
(figure |BJ)- The diagrams are plotted for the same incident wave polarizations as in 
figure El at the frequencies of polarization-split doublet. One can see immediately that 
the transmitted wave polarization is modified considerably if gyrotropy is present. The 
structure still works as a polarizer in a sense that at both peaks the transmitted wave is 
always polarized similarly for any input polarization. This can be seen by comparing the 
graphs in figure EH-c as well as in figure |6ji-f. However, the output polarization states 
are found to depend on the optical activity strength present. This dependence mainly 
manifests itself as rotation, so the output polarization states remain visibly orthogonal 
for the two peaks. However, besides the change in orientation there is also a change in 
ellipticity, as can best be seen in figure HJ. 

It is also worth noting that the character of polarization rotation in this case is 
completely different from what was observed in the previous section. This difference is 
best seen comparing figure|Hfc,f and figure|Hfc,f, when the same incident wave polarization 
is considered, and plotting the angle of rotation versus optical activity a (figured). 
In the bi-isotropic case the rotation is always with respect to incident wave, it is 
independent of the structure's geometrical or spectral features, so <fi is almost equal 
for resonant and off-resonant frequencies, exhibiting only a slight frequency variation 
according to |T0|) (figure [7^). 

In the case when isotropic optical activity is added to birefringence, the 
dominant feature is polarization-induced splitting, its polarization properties subject 
to modification due to gyrotropy. This happens in a multilayer because splitting (along 
with all spectral features) is resonant and the more layers there are in a multilayer, 
the more pronounced the peaks are. On the other hand, gyrotropy-induced rotation 
in (}26|) only weakly increases with the number of layers. So, the rotation occurs with 
respect to eigenpolarizations of the birefringent medium. The angle of rotation is thus 
strongly dependent on frequency, exhibiting a change of 90° over a very small change 
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(Aa;/a;o = 0.0015) between two polarization-split peaks (figure Ujp). The dependence 
4>(a) itself no longer conforms to (jlOj) and even exhibits a slight nonlinearity, apparently 
connected to the emerging ellipticity in output polarization (see figure EF)- 

It may appear at first that this change of output polarization should correspond 
to the change of polarization eigenstates in gyrotropic media. However, in this case 
one would expect exactly the same propagator structure as in (J2IJ) with a and b 
replaced with elliptical eigenvectors, which could have been recoverable from figure El 
This, however, is not straightforward (see |Appendix A\ , and this would also contradict 
with the structure in (|27j) where a term proportional to da is present. Together with 
an additional term in (}26|) this means that eigenvectors should be thickness-dependent. 
This cannot happen since eigenvectors of P are the same as those of M, which depends 
only on material parameters. 

Indeed, if only change of eigenpolarizations were involved here, then gyrotropy 
would leave unchanged the Malus-like dependence of transmittance on input 
polarization, with respect to rotated output polarization states - quite contrary to what 
is seen in figure Eh where the output and input waves are almost mutually perpendicular 
in polarization, yet the transmittance is near unity. So one may guess that there should 
be an interplay between "bi-isotropic-like" uniform polarization rotation and gyrotropy- 
induced change of eigenpolarizations. The exact nature of this interplay is yet to be 
understood. 

To summarize this section, we have found out that isotropic optical activity added 
to a birefringent multilayer structure influences both the transmission spectra and the 
transmitted wave polarization. However, the modification of the spectra can be seen in 
the change of polarization sensitivity. The nature of polarization-induced splitting, the 
frequencies of split components, and the orthogonality of output polarization for these 
components remains unaltered. 

The results obtained can be used to control the output polarization of a splitter 
based on an NPD multilayer with anisotropic layers. This control occurs without any 
modification to the splitter frequency characteristics. This is definitely of interest 
in the design of integrated optical polarizers, as both polarization directions can 
be rotated without having to physically re-orient the multilayer. Again, fractal 
structures are geometrically preferable as NPD structures because peak multiplicity 
allows polychromatic filtering properties. This can facilitate both multiple-channel 
polarization separators and split component matching, which can in turn be promising 
in applications. 

5. Conclusions and outlook 

In conclusion, we have analyzed the influence of isotropic optical activity, the simplest 
kind of gyrotropy, on the spectral and polarization properties of deterministic non- 
periodic (in particular, fractal) multilayer structures. In isotropic case, we have shown 
numerically and analytically that optical activity does not cause any change in the 
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transmission spectrum, but rather leads to a uniform polarization rotation. It is 
independent of the structure geometry and can be described by (|10p. In birefringent case, 
we have shown that the frequencies of polarization-split doublet do not change. However, 
both the output polarization and the value of transmittance at peak frequencies do 
depend on the strength of gyrotropy. In more complicated cases, when gyrotropy- 
related terms are expected to cross-couple with e and /i in the single-layer characteristic 
operator (jlHSJ), more significant modification of the spectra can be expected, e.g., 
changes of resonant peak locations. 

The results obtained are characteristic to any structure geometry that facilitates 
sharp resonance peaks located in a band gap, as is the case with many deterministic non- 
periodic media. Fractal geometry considered here allows the effects to be manifest in the 
desire frequency range with desired peak multiplicity. Among other notable geometries 
one can also name single- and multiple-cavity Fabry-Perot resonators embedded into a 
periodic multilayer. 

The effects observed can be used in the design of frequency-selective, compact 
devices for polarization-sensitive integrated optics. This particularly concerns mono- and 
polychromatic absorptionless polarizers with controlled output polarization. Devices 
that rotate the polarization of output beam with respect to the input are also possible. 
Note, however, that in this paper no account was taken for present-day experimental 
or technological applicability of the results obtained. Also note that the combination 
of parameters described in Section 0] is not among those naturally occurring in crystals 
(though it may be possible in nanocomposites). Nevertheless, this combination is very 
useful illustratively, since it allows analytical treatment and so can be used as a landmark 
of what will happen if we complicate the system further. 

There is a broad range of problems to be addressed within the scope of extending the 
present research. First and foremost, it appears fruitful to investigate more complicated 
cases of optical anisotropy and gyrotropy. This is especially relevant is two cases. First, 
of interest are materials with properties subject to external tuning, e.g., liquid crystals 
or Pockels, Faraday, or Kerr media. Secondly, as was mentioned above, of importance 
are anisotropic and/or gyrotropic materials which can be readily used in multilayer 
micro- and nanostructure fabrication. Since the numerical techniques used in this paper 
are based on the general material equations ([1}, it appears possible to rigorously allow 
for all kinds of anisotropy and gyrotropy, as well as to account for polarization effects 
during wave propagation. This paper demonstrates that analytical or semi-analytical 
treatment is possible using symbolic computation techniques. A thorough classification 
of spectral and polarization-related effects resulting from anisotropy and/or gyrotropy 
according to different crystallographic symmetries existent in nature is worthwhile, too. 

Other areas of extending this paper are also apparent. Namely, it is interesting to 
find out what happens in NPD multilayers containing anisotropic or gyrotropic layers 
of the same material but different spatial orientation. In addition, it also appears 
promising to combine materials with various optical properties in different parts of the 
same NPD multilayer. Since such structures generally tend to exhibit eigenmodes with 
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distinct localization patterns, one may expect that by introducing different materials 
into different localization regions one can produce devices with unusual and largely 
tunable spectral dependencies of optical properties. 
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Appendix A. Bi-isotropic media and circular polarization 



From the reasoning of Section |31 one may conceive that bi-isotropic media allegedly do 
not discriminate between left-handed and right-handed circular polarization, while it is 
commonly known (see [27,33]) that on the contrary, change of effective refractive index 
occurs. Here we would like to explain in detail how and why this effective index change 
does not cause any polarization splitting. 

We take a single bi-isotropic layer and analyze what happens to a circularly 
polarized incident wave. Starting with (fTTj) . we can see that in this case (0 = kda) 



h; 



# (a + ib c 



HnL 



1 



[(cos + i sin 0) a + i (cos + i sin 0) b c 



L 



1 

no 



H V 



(A.l) 



We can see that there is no change of polarization but rather a change of phase for 
the transmitted wave with respect to that for isotropic layer. The value of this change 
is 0. If one reverses the direction of circular polarization, one gets 
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As can be noticed, the change of phase is different in this case. A conventional 
approach with bulk bi-isotropic media is to combine this phase shift with propagation 
phase ip = kdy/ep. This results in a change in effective refractive index, such that [27] 
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rin ± Sn = y/ejl ± a. 



(A.3) 
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However, our studies show that such an approach can yield misleading results when 
applied even to a single layer. A straightforward conclusion of the effective index concept 
would be that the layer, and even more so the multilayer structure, should become 
sensitive to the orientation of circular polarization. Indeed, as spectral features strongly 
depend on n, the difference in n e g should result in the difference in the spectra, causing 
polarization-induced splitting with respect to circular eigenpolarizations. 

However, this contradicts with both numerical and analytical results. Indeed, note 
first that the phase change in (|A.1HA.2|) . even if included in L, does not contribute 
to transmittance. Further, considering a multilayer composed of N isotropic and bi- 
isotropic layers, one can use fT8H20j) to see that 



Since transmission spectrum is determined by the matrix product in (|A.4jl . which 
is polarization-independent, bi-isotropy in any of the constituent layers cannot lead to 
polarization-induced splitting. Instead, it contributes to an overall uniform phase change 
- similar to a uniform polarization rotation in the case of linear input polarization. 

This example is a good confirmation of the fact that the refractive index itself is no 
longer straightforward nor safe to use once complex optical media are involved. Indeed, 
an effective refractive index for a certain eigenpolarization essentially means that a wave 
with this polarization propagates as if the medium were isotropic. But it can only be true 
if the wave polarization is never changed during propagation. In homogeneous media 
of any kind, it is always true by definition of an eigenwave [33] - so it is enough to 
make sure that the incident wave is properly polarized, which is what is conventionally 
assumed. 

However, in inhomogeneous media this is no longer enough. Even in a multilayer 
the layer interfaces can and do introduce polarization coupling between eigenwaves. And 
this is exactly what occurs in a bi-isotropic layer - both layer interfaces serve as perfect 
couplers between eigenpolarizations, reversing the polarization after each reflection. 
(This cross- coupling is mentioned in [26] and is due to the fact that right- or left- 
handedness of a circularly-polarized wave is determined by the sign of the vector product 
in[H x H*] [22], and during Fresnel reflection in our case n — ► — n, H — > rH, r being 
real). As such, even in the case of one layer, even in the case of normal incidence and 
axial symmetry of the system, the wave initially polarized along one of the eigenwaves 
does not retain its polarization during propagation. Therefore the resulting spectral 
properties even of one layer differs greatly from what would have happened if one blindly 
used the effective index concept suitable for homogeneous media. 
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Appendix B. Evolution operator for a uniaxial, isotropically optically active 
layer 

Here we would like to show in detail the way we have derived the model propagator for 
a uniaxially birefringent layer with isotropic gyrotropy to be used in (J2SJ) and (|2*7j) . 

The main difficulty in evaluating the propagator directly from f!23[) consists in the 
fact that eigenvectors of such a matrix are complex elliptical vectors whose orientation 
and ellipticity depend on numerical relations between e , e e , and a, and layer interfaces 
even for one layer will play a role not clearly understood (see |Appendix A| ). The 
expression for P becomes excessively complex to be analyzed successfully. 

One of the ways out would be to construct approximations assuming that a and 
|e e — e G | are small quantities. However, such approximation has to be done carefully, 
since while these quantities are small enough indeed, they appear in the propagator 
multiplied by du/c, and the resulting phases are no longer infinitesimal. So it can be 
misleading to straightforwardly suppose \e e — e \ 2 = = a 2 . 

Instead, what can be done is to analyze several orders of Taylor series for P in order 
to recover how the Taylor decomposition can fold back into trigonometric expressions so 
as to produce the same Taylor series up to a certain order. The result of this operation 
is what we call here the model propagator. Reduced to a 4 x 4 form [see the footnote 
for (Oil , it reads 
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In order to simplify off-diagonal terms to the maximum possible extent, we have 
substituted e Q = e — Se, e e = e + 5e. Also, A = | sin(<5e/e) and e' = e (1 — 8<p) where 
5 tp = 5e(kd^/€jl)~ 2 , and as before, (p = kda. 

It can be confirmed that the propagator constructed according to (jB.lHB.5|) yields 
the same Taylor decomposition as the real matrix exponential when the terms up to a 2 
and Se 2 are preserved. 

Certainly, the propagator of (|B.lHB.5]l is by no means a true form of the evolution 
operator; it is not even an approximation in the strict sense - rather, it is an expression, 
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which yields the same approximation. Hence it can be used as an approximate model 
of a real gyrotropic layer. The fact that it is trigonometrical allows us to analyze such 
phenomena as polarization plane rotation and birefringence, reverting to decomposition 
only when and where needed. The only apparent limitation of this approach is that 
we should not decompose (jB.2UB.5j) to more than the second order with respect to a 
and Se. 

This kept in mind, one can take Taylor in square roots with respect to Se and neglect 
the product aSe (which is essentially a quadratic term). This yields (ip = kd^/ep) 

A = D ps P cos ip + SeXdS sin ip, 

B ps -f = [P sin if - Se (X od + A') S cos cp) , 

C ps [P sin ip + Se (\ od - A') S cos <p\ , 

P = I cos (p — q x sin (ft, S = a (8> a — b <8> b 

where the subscripts d and od stand for "diagonal" and "off-diagonal", and 

kdji kd^lj 1 

A = -. (B.7) 

analogous to (jTlj) and (J22J: 

(B.8) 



;b.6) 



A n 



A 



Eq. (JR 



2^e/i 2e 
Therefore one can rewrite (jB.6MB.7jl in the form analogous to (jTlj) and (j2H 
P (9) ps L (/ cos - q x sin 0) + <5eL' (a ® a - b ® b ) 
is used in the text as (l26l). The coefficient matrices are 
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On the other hand, one can take Taylor of (jB.2HB.5j) in another manner, 
excluding the propagation phase ip rather than the polarization plane rotation <fi from 
decomposition. Making use of the fact that 

ie . , , i — ; — ; — rr % (e ± Se) 



sin kd\/ (jl (e ± Se) 



(e ± Se) 



Vt 1 ( e =F Se ) 
and likewise neglecting aSe, we obtain 
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structure on analogy with (|B.8[) . we have 

P (9) ps L o a <g> a + L e b ® b - /ccfaLq* . (B. 12) 

This formula is used in the text as (j27j) . The coefficient matrices L a and L& are 
identical to those in (j25j). The antisymmetric term contains the same matrix L as 
in (E 
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